First-order coherent resonant tunneling through an interacting coupled-quantum-dot 
interferometer: generic quantum rate equations and current noise 
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We carry out a detailed analysis of coherent resonant tunneling through two coupled quantum 
dots (CQD) in a parallel arrangement in the weak tunneling limit. We establish a set of quantum 
rate equations (QREs) in terms of the eigenstate-representation by means of a generic quantum 
Langevin equation approach, which is valid for arbitrary bias-voltage, temperature, and interdot 
hopping strength. Based on linear-response theory, we further derive the current and frequency- 
independent shot noise formulae. Our results reveal that a previously used formula for evaluating 
Schottky-type noise of a "classical" single-electron transistor is a direct result of linear-response 
theory, and it remains applicable for small quantum devices with internal coupling. Our numerical 
calculations show some interesting transport features (i) for a series-CQD: the appearance of a NDC 
due to the bias- voltage-induced shifting of bare levels or a finite interdot Coulomb repulsion, and (ii) 
for a parallel CQD in strong interdot Coulomb repulsion regime: finite-bias-induced AB oscillations 
of current, and magnetic-fiux-controUable negative differential conductance and a huge Fano factor. 

PACS numbers: 72.10.Bg, 73.63.Kv, 73.23.Hk, 72.70.-|-m 
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I. INTRODUCTION 

The investigation of quantum oscillations between two 
levels in a coupled quantum dot (CQD) system under 
transport conditions has been the subject of enormous in- 
terest over the last decadeiii^ Recently, a CQD arranged 
in parallel between source and drain has been experimen- 
tally reported to constitute a mesoscopic QD Aharonov- 
Bohm (AB) interferometerj^i^ and it is believed that 
manipulation of each of the QDs separately and the ap- 
plication of magnetic-flux piercing the device can provide 
controllable parameters for the design of transport prop- 
erties by tuning quantum oscillationsi^'^'^'^'^'^ '^^i^^'"'^^'^^ 

To describe such quantum oscillations in quantum 
transport, master equations and a "quantum" version 
of rate equations were first proposed by Nazarov and 
coworkers jiiiS and later derived from the Schrodinger 
equation,-A§!AZ respectively. These original works are 
mainly for a CQD in a series arrangement between leads 
in the limit of zero temperature and large bias-voltage. 
The authors have generalized the quantum rate equa- 
tions (QREs) for the case of a CQD interferometer in 
the same limit and employed them to study magnetic- 
flux-controUcd photon-assisted tunnelingr^- Recently, a 
bias-voltage- and temperature-dependent generalization 
of the QREs was carried out by the authors employing 
the nonequilibrium Green's function in the limit of weak 
dot-dot coupling. However, the dot-dot coupling of the 
CQD device is usually tuned using gate voltage in ex- 
periments, and it may not be weaker than the tunnel- 
coupling to external electrodes. Therefore, it is desirable 
to develop generic QREs without assuming weakness of 
the dot-dot couphng for the purpose of systematically 
analyzing the transport properties of a CQD system. 



To accomplish this, we employ a generic quantum 
Langevin equation approach to derive the QREs in the 
eigenstate representation. The main result is the deter- 
mination of a new coherent transfer term that emerges 
in all dynamic equations of the reduced density-matrix 
(RDM) elements stemming from the effective coupling 
between two eigenstates due to tunneling, which is absent 
from our previous results for QREs.^* In addition, other 
coherent terms occur in the ensuing QREs for a CQD 
interferometer due to interference between the two path 
branches enclosing magnetic-flux, which is responsible for 
the AB oscillation feature of transport, as expected. 

In another emerging aspect of CQD systems, stud- 
ies of current fluctuations have become an important 
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Several analyses of shot noise in a CQD 



system have been undertaken by means of the QRE 
approach,2ii22i2^i2ii2^i^ZLi2^i^ Most of these works have 
involved calculations in the large bias-voltage limit for 
CQDs in series^i2^iS^i^'2£ or in paralleli^ Our earlier 
work was carried out with analyses based on our previ- 
ous QREs, and, correspondingly it is only valid for the 
case of weak dot-dot hopping. Moreover, some other 
recent studies of bias-voltage-dependent shot noise of a 
CQD have not treated the quantum coherence effect j ^'^i^^ 
Therefore, in this paper, we will also analyze zero- 
frequency shot noise of an interacting CQD using the 
presently developed QREs, focusing our attention on the 
coherence and interference effects and its magnetic-flux 
dependence at both large and small bias- voltages. 

The outline of the paper is as follows. In Sec. II, we 
present our physical model for a CQD interferometer con- 
nected to two electrodes and derive a set of QREs us- 
ing a generic quantum Langevin equation approach with 
a Markovian approximation in describing the dynamic 
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evolution of the RDM elements under transport condi- 
tions. Such microscopically derived QREs are valid at 
arbitrary bias- voltage and temperature, and dot-dot hop- 
ping strength as well. In this section, we also derive the 
current and frequency-independent shot noise formulae in 
terms of the RDM elements using linear-response theory. 
Employing the obtained formulae, we then investigate 
the transport properties of a CQD in series in Sec. Ill, 
stressing the asymmetric transport property. Moreover, 
we study in detail the combined effect of the additional 
coherent transfer term and the two-pathway-interference 
term on current (Sec. IV) and on shot noise (Sec. V) of 
a CQD interferometer. A summary is given in Sec. VI. 

II. HAMILTONIAN AND FORMULATION 

The parallel-coupled interacting QD interferometer 
connected to two normal leads is illustrated in Fig. 1. 
The Hamiltonian is given by 

H = Hl + Hr + Hd + Ht, (1) 

where H^i (rj = L, R) describes noninteracting electron 
baths in the left, right leads, respectively: 

= XI (^"k - /"r,)aj,ka^k, (2) 
k 

with a^k being the annihilation operator of an electron 
with momentum k, and energy £j,k in lead 77, and /i^ 
being the effective Fermi energy of lead rj. In our studies, 
a bias-voltage V is taken to be applied symmetrically 
between the two electrodes, i.e., /xl = —fiR = eV/2, to 
sustain a persistent electron flow from one lead to the 
other. In equilibrium, V = 0, we set 11 l = P^r — 0. 
Throughout, we will use units with h = ks = e — 1. 




FIG. 1; Schematic diagram for coherent resonant tunneling 
through a parallel-coupled double quantum dot system in an 
Aharonov-Bohm interferometer. 

The Hamiltonian of the isolated CQD system, Hdi is: 

Hd= ^ £jc]cj + Uc\cic\c2 + ^{c\c2 + 4ci), (3) 

where Cj is the annihilation operator for a spinless elec- 
tron in the jth QD {j = 1,2). Sj is the energy of the 



single level in the jth QD, measured from the Fermi en- 
ergy of the two electrodes at equilibrium, £1(2) = Sd^^^ 
with 5 being the bare mismatch between the two bare 
levels. Here, we assume that only one single electron 
level in each dot contributes to current. It should be 
noted that since we consider first-order resonant tunnel- 
ing in this paper, and take no account of spin-flip scat- 
tering and spin-flip cotunneling processes (second-order 
tunneling processes), it is reasonable to assume spinless 
electrons in our model. The second term represents the 
interdot Coulomb interaction U. The last term of Eq. ([3]) 
denotes hopping, Q, between the two QDs. 

In this CQD system, there are a total of four possible 
states for the present system: (1) the whole system is 
empty, |0) = |0)i|0)2, and its energy is zero; (2) the first 
QD is singly occupied by an electron, |1) = |l)i|0)2, and 
its energy is £1; (3) the second QD is singly occupied, 
1 2) = 1 0)1 1 1)2, and its energy is £2; (4) both dots are oc- 
cupied, \d) = |l)i|l)2, and its energy is Si + 62 + U. 
Furthermore, with the four possible single electronic 
states considered as the basis, the density-matrix ele- 
ments may be expressed as poo — |0)(0|, pu — |1)(1|, 
P22 = |2)(2|, pdrf = \d){d\, and pi2 = |1)(2|. The statis- 
tical expectation values of the diagonal elements of the 
density-matrix, poo (poo), Pn = (Pn) [j = 1,2), and 
Pdd — (pdd) , give the occupation probabilities of the elec- 
tronic levels in the system being empty, or singly occu- 
pied in the jth QD by an electron, or doubly-occupied by 
electrons, respectively. The off-diagonal term pi2 = (P12) 
describes coherent superposition involving two electronic 
occupation states, |l)i|0)2 and |0)i|l)2. 

To properly account for interference effects between 
the two pathways for tunneling through the system at 
hand, it is convenient to diagonalize the Hamiltonian of 
the isolated CQD by a unitary transformation 

\a) = cos^|l)+sin^|2), (4a) 

1/5) = sin-|l) -cos-|2), 61 = arctan-. (4b) 
2 2 

The transformed Hamiltonian is 

Hd = A„|a)(a| + A^|/3)(/?| + (ei +£2 + U)\d){d\, (5a) 
where Xaip) is the eigen-energy, 

K{i3) = £d ± A, (5b) 



with A = \/Ti? + S^. Correspondingly, the density- 
matrix elements in the new double-dot eigenstate basis, 
Pxx = lx)(xl and pa0 = \a){(3\, have similar physical 
meanings to those in the site- representation (SR). The re- 
lations between these density-matrix elements of singly- 
occupied states in different bases can easily be deduced 
from the unitary transformation, Eq. (jH). (Note that the 
empty state, poo, and the double occupation state, pdd, 
have the same meaning in the both representations.) 



The tunnel-coupling between the interferometer and 
the electrodes, Hx, can be written in the eigenstate- 
representation (ER) as {sg — sin | and cg — cos |) 

x(ce|0)(a| + sg\0){P\ - Sg\a){d\ + Cg\f3){d\) 

x{sg\0){a\-cg\0){(3\+ce\a){d\+sg\(3){d\)+Ii.c.l6) 

For simplicity, the lead-dot tunneling matrix elements, 
Vrij, are assumed to be real and independent of energy. 
The factor e^^'^^^ is the accumulated Peierls phase due 
to the magnetic-flux $ ((^ = 27r$/$Oi ^'o = hc/e is the 
magnetic-flux quantum) which penetrates the area en- 
closed by two tunneling pathways of the interferometer. 

In the following, we apply a generic quantum Langevin 
equation approacb22i^i^i2^i^i2^i^ to derive a set of 
quantum rate equations (QREs) to describe the dynam- 
ics of the system variables of the CQD due to coherent 
resonant tunneling between external reservoirs, as mod- 
eled by Eq. H]). In the derivation, three steps are in- 
volved: First, we start from the Heisenberg equations 
of motion (EOMs) for the density-matrix operators pooj 
p^^i , and pdd in the ER and related reservoir operators 
c,,k, and then formally solve them by integration of these 
EOMs exactly. Next, under the assumption that the time 
scale of decay processes is much slower than that of free 
evolution, which is reasonable in the weak-tunneling ap- 
proximation, we replace the time-dependent operators in- 
volved in the integrals of these EOMs approximately in 
terms of their free evolutions. Thirdly, these EOMs are 
expanded in powers of the tunnel-coupling matrix ele- 
ment Vjjj up to second order. By adopting a Markovian 
approximation, we finally derive the generic QREs with 
arbitrary bias-voltage and temperature, as well as arbi- 
trary dot-dot hopping, as 

POO = -[c^(r^i„ + r+^) + s2(r+^^ + r+j 
+- sin6i(r+„ + r+„ - - r+^)]poo 

+ [ce^iia + 4^22a + 2 sin^(rr2a + ^21 J]P»o, 
-f- sin6'(rj;i„ + r-^^ - T:^^^ - T-^p){pafj + pi3a) 

+ 2^4{^12a + ^I2p) - cl{T21a + r^l/s)]^"/? 
-7^[cl{^12a + ^I2f}) - 4{^21a + r2i^)]P/3a, (7a) 



Paa = [cer7"i„-i-s0r+„-t--sin6'(r+„ + r+„)]poo 

-[cl^lla + 4^22a + \ siu 6'(rj;2„ + r2"i„)]paa 
~[4^11I3 + '^1^2213 ^ 2 ^^'^^(12/3 + ^2ll3)\Paa 



'2^4(^120 ^ ^12a) ^ 4(^21fJ ^ ^21a)]Pa0 
'i[4{^2ip - ^tla) ~ cli^Ufi ~ f +^)]p;5„, (7b) 



hp = [^er+^ + CerJ,/?- ^sin0(r+^-|-r+^)]poo 

~\4)^\i(i + 4^22Si - 2 sin6'(rj;2^j + ^2i!i)\PPSi 

-[c^f+„ + s2f+^ + i Sin0(f + „ -f f^+i J]p/50 

-- sin6'(r~ ^ - r22„ - r+ ^ + r+ ^)(pa^ -i- p^^) 

-\\4^21a - fJl/j) - Ce(r-2^ - f +^)]p/3a, (7C) 



Pdd — \4^X\jj ^ 4^X2^ 2 ^^'^'•^12/3 ^ "'^21/3^]^"" 

+ [c^f + „ + S^f +2 „ + \ Sin(f +2„^+ f + J]p3/3 

"[■59(^11^ + ^22^) + 4^\\cx + ^22/3) 

+ 2 '5in6'(r72„ - ^^2(i + ^2\a - ^2\ti)\Pdd 

sin6'(r+ „ + - r+ „ - rJj^)(pa/3 + Pfia) 
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¥An2a + r+fl) - 5^(rJi„ + r2+ «)]p,3„, (7d) 



= f2Ap„^ -(- [- sin f?(r+ „ + r+ - r+2„ - r+ ) 



2c^(ri2a + r^2;3) + T^^^^tia + r2+i^)]poo 



2 V^e^Wa + -1^22/3 + ^X\a + ^X,2(i) 
^4^\\(} + ^220 + ^X\(i + r^2a) 
+ 2 ™^(rr2a + ^2\o, - ^120 ~ ^210 

+^120 + tJlq ^ rj*2^ — rJ^^)]pQ/3 

+ 2['^e(^12a ~ ^X2(i) ~ '^9(^210 ~ '^X\(i) 
-i sine'CFj^i^ - r22a - f + f j2/3)]Paa 
+ 2 [^^6(^12/3 ^ ^X2a) ^ *e(r21/3 ^ rj;^^) 

-i sin6'(r7^^ - r22^ - f + „ + r+2„)]p,3,3 

+ ^[^ sin6'(f 22„ + r22^ - f j^i^ - P-^^) 

+Ce(rr2a + "^1213) - Se(r21a + f 2"i^)]Pdrf, 



(7e) 
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where the electron tunnehng-in(out) rates are defined as 



i 12x 



Er^i^x = Er.i2e^"^'^/V,f (Ax), (8b) 



33X 



21x 



Ef.1« = Er.../.^(^x + c^)- 



(8c) 
(8d) 



Er.^i2x = Er.i2e^"^'^/V,f (A, + C/), (8e) 

ri ri 

Er^2ix = Er.i2e^^''^'^/'/,^(Ax + u), 



(8f) 

[/,j'(e) is the Fermi distribution function of lead rj, 
f-{e) = 1 - /+(e), and s^/;? = ±1] with 



(9) 



being constant in the wide band limit (p^ is the density of 
states of lead 77). T^jj' represents the strength of tunnel- 
coupling between the jth QD level and lead 77 if j = j', 
and otherwise (j ^ j') it measures the interference in 
tunneling through the different pathways. The adjoint 
equation of Eq. (fZel) gives the equation of motion for the 
RDM off-diagonal element piSa- Also, we must note the 
normalization relation poo + Pxx + Pdd = 1- Simple 
rate equations written in the ER were obtained earlier in 
Ref. [33 for double-well semiconductor heterostructures. 

In the ER, the CQD interferometer is equivalent to a 
single QD with two energy-levels, both of which couple to 
the two electrodes such that either level is tunnel-coupled 
to either electrode via two different pathways with differ- 
ent tunneling matrix elements depending on differing ac- 
cumulated phases due to the magnetic-flux [see Eq. 
Figure [3] below exhibits schematic diagrams of energy 
configurations for coherent tunneling through a CQD in 
the ER. In this situation, all tunneling events can be 
classified in three different categories. For example, F^;^^ 
(r||"-^^) describes the tunnefing rate of an electron enter- 
ing (leaving) level a in the CQD without the occupancy 
of level (3 via pathway 1 (indicated by the upper line in 
Fig. 1); Similarly, F^aa (^220) denotes the same tunnel- 
ing process via pathway 2 (indicated by the lower line in 
Fig. 1); T'^2a (^i2a) dcscribcs the interferential tunneling 
process of an electron entering (leaving) level a without 
the occupancy of level /? due to the interference between 
the two pathways 1 and 2. All ^fj,^ are the correspond- 
ing rates in the case of double occupation. Naturally, 
only the interferential tunneling term suffers an accumu- 
lated phase factor as shown in Eqs. ([5c)) and ((8^ . 

Therefore, the classical parts of the dynamical equa- 
tions of the RDM diagonal elements have clear classical 



interpretations. For example, the rate of change of elec- 
tron number in level a, Paai in the CQD, governed by 
Eq. (j7bp . is contributed to by the following four single- 
particle tunneling processes: (1) tunneling into level a of 
the CQD from both left and right leads, if the CQD is ini- 
tially in the empty state poo; (2) tunneling out from(into) 
level a{l3) of the CQD into(from) both leads, when the 
CQD is initially just in the state Paa{pi3p)] (3) tunnel- 
ing out from level /? of the CQD into both leads, when 
the CQD is initially in the state pdd, via path 1, path 
2 and interference contributions. The last three terms 
in Eq. (j7b[) describe transitions between two eigenstates 
through the effective coupling with off-diagonal elements 
via both paths and interference, which have no classical 
counterpart. They are responsible for coherent (quan- 
tum) effects in the transport. We note that even though 
there is no direct coupling between the two eigenstates 
of the isolated Hamiltonian, Eq. ([5a|) . in the ER, the 
tunnel-coupling described by Eq. ^ results in an effec- 
tive transition between them and thus leads to a quantum 
superposition state between the two eigenstates, whose 
dynamics is ruled by the equation of motion of the off- 
diagonal element pc^, Eq. (|7e)) . 

The tunneling current operator through the interfer- 
ometer is defined as the time rate of change of the charge 
density, = Eka!,ka,,k, in lead 77: 



(10) 



According to linear-response theory in the interaction 
picture, we have:^^ 



Il = {Jl) = -i dt'{[JLit),Hjit')U. (11) 

J —00 

Along the same procedure as above, we arrive at the fol- 
lowing expression for the current within the weak tunnel- 
coupling approximation, 



+ - sin6'(Fj^2Q + ^L21a ~ ^L12p ~ ^L21l3)] 



-[('l^Lllc + 4^L22c. + 2 '5in6'(F^^2^ + F^2^^)]p„„ 
'[^l^llip + '^I^L22P ~ 2 ^i'^^(-'"L12;3 + '^L21f})\Paa 
^[4^lllfS + C^9^l22fi - \ sin6'(rZi2/3 + ^l2Ui)\Pl3P 

~2^2 ^^'^^^"'^^11" ^ ~ -'^L22a ^ ^ L2213) 



'^l(^L21c 



- L21l3)\PaP 



- F+ - F+ 1 
- L22/3 ^ LllQ ^ L11I3> 



+ *e(rL12a + ^ L12p> 

+ i[lsin0(f+... +F^ 

"■Se(rL12a + ^ll2fj) + Ce(rj2lQ + ^i2ll3)\PaP 
~2^2 ^"^"^^^^Llla + ^Lllf} ~ ^L22a ~ ^ L2213) 
~^e(^L12a + ^^12/3) + Se(r^21a + ^ L21l3)]PPa 
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+ Ce(rLi2Q, + ^tl2fi) ~ *e(£L21a + ^L2ip)]PPa 
~{^e^Llip + ^e^Llla + *9rL22a + '^0rL22/3 

-- sin6'(r^^2/3 - r^^2a + ^L21f3 ^ r^2lQ)]Pdd- 

(12) 

Interchanging the subscripts "L" and "R" , we obtain the 
tunnehng current In relevant to the right lead. It is easy 
to verify that II = —Ir- This current formula demon- 
strates that all possible tunneling events relevant to lead 
r] (tunneling through both paths and with the interfer- 
ential term) provide corresponding contributions to the 
current of lead 77; the current is determined not only by 
the RDM diagonal elements, but it also involves the off- 
diagonal elements explicitly. 

Another observable of interest is the current noise, 
whose spectrum is defined as the Fourier transform of 
the symmetric current-current correlation function 

r°° 1 

5wH = y dTe^'^^-{[SJr,{t),SJ,,it')]+), (13) 

with 6Jri{t) — Jri{t) — Irj- Here, we only consider the 
frequency-independent (Schottky-type) shot noise. Em- 
ploying linear-response theory we can also evaluate its 
statistical expectation value to the first nonvanishing 
term in iJ/, obtaining 

Sll = \ j°^jT{[Jl{t),Jl{t')]+) 

= ['^l(^llla + ^^22/3) + ^li^tllfi + ^^220) 

+ 2 sin ^(rJi2a + rj21a " ^tl2p ' ^^21/3)1^00 

+ [(^i'^lllcc + SerZ22a + 2 Sin6l(r^i2a + ^l2la)]Po.o. 

+ [^1^11113 + '^0^^22/3 ~ 2 '^^'^^(^il2/3 + ^L2w)\Paa 

+ [s6^lll(3 + CerZ22/3 - \ Sinf (rZi2/3 + rZ21/3)]P/3/3 

+ 2[2 +^ill/3 ~ r^22Q ~ ^L22li) 

+^e{^Ll2a + ^L12I3) ^ ^e{^ L21a + ^ L2lf3)\PaP 
+ Sin0(fj22a + ^122P~ "^lll. ' ^lll^) 

~^l{^ll2a + ^Il2f3) + '^e(rL21a + ^l2lfi)\Paf3 

+ 2^2 ^^^■^^^^^11" ^ '^Llip ~ -'^L22a ~ ^L22p) 
~Ce(r^l2a + '^L12p) + ^l)i^L21a + L21l3)\Plia 

sin0(r+22„ + r+22^ - r+ii„ - r+^i^) 

+^l(^ll2a + ^^12/3) ~ ■*e(£L21a + ^t2lf3)\PPa 

'^[^e^Liip + "^er^j^Q, -I- SQTj^22a + '^er^22;3 

-- sin6l(r-^2,3 - f Zi2„ + r£2i0 - ^'L2lJ]Pdd- (14) 



The corresponding result for Srh is obtained by inter- 
changing the subscripts "L" and "R" in Eq. (fT4| . In 
leading order, Slr = Srl = 0. 

This shot noise formula merits further discussion. It is 
well-known that the Schottky noise originates from the 
self-correlation of a given tunneling event with itself. In 
literature, this term has been elaborately studied for a se- 
quential picture of resonant tunneling through a double- 
barrier transistor using a master equation^i^ and also 
using a stochastic wave function approach. These pre- 
vious studies show that the "classical" intrinsic shot noise 
can be simply expressed by the formul a^^i^° 

n 

provided the current can be written in the following form 
^'; = E(+/-)Wn: (15b) 



where p is a vector whose components are the RDM el- 
ements of the single-electron transistor, the superscript 
"0" indicates its stationary solution based on the rate 
equation, and 7^„ is the corresponding rate of a possi- 
ble tunneling process that can change the state n of the 
single-electron transistor and contribute to the current 
of lead 77. This is to say that all single-particle tunnel- 
ing processes involving lead 77 will give a positive (-I-) 
or a negative (— ) contribution to the current of lead 
rj, depending on whether an electron is to enter or to 
leave this lead in a particular tunneling process; never- 
theless self-correlation terms in all these same processes 
will always result in a positive (-I-) contribution to the 
Schottky-type shot noise, and the only difference be- 
tween the current and the self-correlation noise formulae, 
Eqs. (jl5bp and (|15ap . is a sign change. We note that the 
formula Eq. ()15ap has been only proved for a "classical" 
double-barrier tunneling device, in which the RDM of 
the device has only diagonal elements (thus p° > 0) and 
the tunneling rate jrjn is always non- negative (7,,„ > 0). 
Without attempting to rigorously verify its validity, this 
scheme has been recently applied directly to calculate the 
Schottky-type shot noise of a QD system having coupled 
internal degrees of freedom, but with its current being ex- 
plicitly dependent only on RDM diagonal elements which 
always have positive values ] ^^'^^ The present theoretical 
formulation circumvents this problem. Furthermore, we 
prove here that (1) this scheme is just a result of linear- 
response theory in the weak-tunneling limit; and (2) this 
scheme remains valid for the calculation of Schottky-type 
shot noise of a CQD interferometer, for which not only 
RDM diagonal elements contribute to the current, but 
also off-diagonal elements contribute as well, with sta- 
tionary values and corresponding tunneling-rates which 
could possibly both be negative. 
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III. CQD IN A SERIES CONFIGURATION 

In this section, we discuss coherent resonant tunnel- 
ing through a series-CQD in the case of infinite inter-dot 
Coulomb repulsion, U = oo. In this case, the correspond- 
ing tunneling amplitudes and the interferential tunneling 
vanish, Tl22 = '^mi = 0, T^ia = Tri2 = 0, and all 
r = 0. We also choose 6 — {ei = e2 = Ed)- 

For comparison with our previous QREs in the case 
of weak inter-dot hopping,-i^ we re-express the resulting 
generic QREs Eqs ([7a| - ([7e|) in the SR for the case of a 
series configuration (rf = rf^^, and r| = r|22^): 



Poo 



Pii 



rt)poo 



^lPoo ~ r^/Oii + if7(p2i - P12) 



P22 = r+poo - r^P22 - i^{p2\ - P12) 

--sin6i(r^„ -r^^)(pi2 +P21) 



P21), 

(16a) 
(16b) 

(16c) 



as shown by the final three terms on the rhs of Eq. (|16dp . 
Moreover, the real part of the off-diagonal elements make 
an explicit, direct contribution to the tunneling current 
and Schottky-type shot noise as well. 

In regard to the fact that our previous QREs are only 
valid for the case of weak interdot hopping, whereas our 
present QREs are derived without that limitation, it is 
interesting to note that the present generic QREs do re- 
duce to the previous QREs in the limit — > (^ ^l/r)- 
In this situation, from Eq. (|5b[) we have Xa(ff) = i 



and thus 



This 



(S = 0), leading to Aa « A/5 c,,^ -xx^^o ^ ~ x,^^. 
approximate relation eliminates the additional quantum 
coherence terms of Eqs. (|16ap - ()16dl) due to the tunneling- 
induced effective dot-dot hopping, and consequently the 
QREs reduce exactly to our previous results. For numer- 
ical verification of this assertion, we examine the effect 
of hopping (ri/r — 0.2 1.0) on resonant tunneling in 
Fig. [21 where we plot the measurable physical quanti- 
ties, e.ff. the occunation numbers in both dots, n^^ a,nd 



P12 = i2Spi2 - i^ipii - P22) - -^(Tl +^R)Pi2 



+ ^sin0(r+„ 
-7sin0(r^„ 



r+ 

^ Ra 



Lf3 



Rf3 



I Poo 



T 



RpjPii 



4sin0(r^, 



r 



with Poo + Pii + P22 = 1 and p2i = P12; in which 



and 



p± _ 2p± 
L — '^0^ La 



p± _ 2p± 

R — Ra 



Lp)P-22' 

(16d) 



(16e) 



(16f) 



denotes the effective tunneling rates between the 
first(second) dot and the left(right) lead. The current 
becomes 

II = r+poo-rzpii-^sin6i(r^^-r^^)(pi2-fp2i). (i7) 

Comparing with our previous results [Eqs. (39a)- (39c) 
of Ref. [1^, we find some new features in our present 
QREs, Eqs. lfTCa|) - (pd| . in the SR: (1) The tunneling 
rates between dot and lead, F^^^, are modified to ac- 
count for the effect of interdot hopping; (2) More in- 
terestingly, effective interdot hopping induced by tun- 
neling to external electrodes provides a new contribu- 
tion, which is dependent on the real part of the RDM 
off-diagonal elements, to the dynamics of the RDM di- 
agonal elements stemming from quantum coherence [for 
instance, the fourth terms on the right hand side (rhs) of 
Eqs. (|16b|) and (|16cp ]: (3) Conversely, this effective inter- 
dot hopping also leads to a new damping of quantum su- 
perposition related to the diagonal elements themselves, 
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FIG. 2: Average current (a,c), and the differential conduc- 
tance (b,d) vs. bias-voltage with several different interdot 
flopping strengths Q,/V = 0.2, 0.5, and 1.0. (a) and (b) are 
plotted for Ed/T = 1.0, (c) and (d) for e^/r = -1.0. The 
thick fines denote the resufts cafcufated by the present QREs, 
Eqs. (|16a|l - (ll6d[) and Eq. HZ]); whife the thin fines arise from 
the previous QREs, Eqs. (39a)- (39c) in Ref.[ii. The temper- 
ature is fixed as T/F = 0.2. 

It is obvious that there is no nontrivial difference in the 
separate calculations for these physical quantities in the 
case of weakest dot-dot hopping, fi/F = 0.2. With in- 
creasing dot-dot hopping, the difference between the two 
theories becomes unambiguous and cannot be ignored. In 
particular, a zero-bias peak (ZBP) in the differential con- 
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ductance, dV/dl, is found by the present theory for the 
system of Sd/T — 1.0 with dot-dot hopping Q/T = 1.0, 
in contrast to the peak-sphtting around zero-bias- voltage 
predicted by the previous QREs. On the other hand, for 
the system of e^/r — —1.0, both calculated results for 
differential conductance show a peak-splitting structure 
but they have different splitting widths. The transport 
properties are clearly asymmetric between the systems 
with Ed > and Ed < 0, which have been already pointed 
out in linear transport regimei^ In fact, this asymmet- 
ric feature can be easily explained by examining the en- 
ergetic structure of the present system. As mentioned 
above, the CQD in the ER can be mapped onto a model 
of a single QD with two levels, a and /3, each connecting 
to both electrodes with different tunneling matrix ele- 
ments. The physical picture is illustrated in Fig. [31 which 
includes a total of six configurations [(a)- (f)], showing the 
relative relations among the two eigen-energies and the 
Fermi energies of the two electrodes. For simplicity, we 
set T = in the following analysis. In configurations 
(a) and (d), the CQD is always singly occupied by an 
electron, while it is always empty in configuration (b). 
The three cases have all vanishing first-order tunneling 
current and current noise because the two eigen-levels are 
both far away from resonance with the Fermi levels of the 
two electrodes. On the contrary, both eigen-levels are in 
resonance with the Fermi levels of the electrodes in con- 
figuration (c), which is actually equivalent to a picture of 
resonant tunneling under extremely large bias- voltage. It 
is therefore not unexpected that the resulting rates and 
current in this case are identical to the results obtained 
by Gurvitz and Prager for the same system,-lS as well as 
matching our previous results for the same conditions.^® 

Moreover, it is quite surprising that the current is also 
zero in the resonant configuration (e). We can ascribe 
this result to the strong Coulomb blockade effect, which 
blocks the entry of an additional electron into the singly- 
occupied CQD even though the upper eigen-level is in res- 
onance with the Fermi levels of electrodes. The situation 
is different in configuration (f), in which only the lower 
eigen-level is in resonance with the Fermi levels of two 
leads and thus a nonzero current occurs. We note that 
this current is always smaller than the current in config- 
uration (c) [This is only true for a series-CQD. While for 
a parallel-CQD, the situation may be opposite (see be- 
low)]. Actually, configuration (f) describes a picture of 
resonant tunneling under a small bias- voltage, in contrast 
to the case of large bias- voltage in configuration (c) . 

It should be noted that our present QREs are also 
applicable for the cases with a detuning of the bare 
level energies, i.e., S Q. Actually, spatial variation 
of bias-voltage on the series-CQD is quite a central is- 
sue in the present studies4^i^i^ Without loss of gener- 
ality, we assume that the bare QD energies are shifted 
by the bias-voltage as £1(2) = ± ijV, and assume fur- 
ther that one fourth of the applied voltage drops at the 
contacts between the QDs and the leads, i.e., the bias- 
voltage- induced shifting factor r] — 1/4. As an illustra- 
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FIG. 3: Schematic energetic configurations for coherent tun- 
neling through a CQD in the ER. The CQD behaves as a 
single QD with two levels a and /3, both of which are coupled 
to both the left and right leads. In configurations (a)-(c), the 
additional quantum coherence terms automatically vanish in 
the QREs (|16a|) - H16d|) and in the current formula (|17p . while 
in (d)-(f), they play an important role in tunneling. 



T/r=0.05 



E/r=i.o 
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n/r=i/2 




0.4 r 



(b) 




FIG. 4: (a) Occupation numbers pn (solid lines), P22 (dashed 
lines) and (b) average current vs. bias-voltage with the bias- 
voltage- induced shifting factor i] = (thin lines) and 1/4 
(thick lines) for ea/T = 1.0, Q/T = 1/2, and T/T = 0.05. 



tion, we exhibit the effects of energy shifting on the elec- 
tron occupation numbers of two QDs and the current 
in Fig. m It is found that the main effect of the spa- 
tially varying electric field is to cause an appearance of 
negative differential conductance (NDC) in the current- 
voltage characteristics, which has been reported previ- 
ously by nonequilibrium Green's function calculations^ 
in order to explain the experimental results in the Tour 
molecules. Here, we can give a simple interpretation 
for the appearance of NDC as follows: Due to the spatial 
variation of bias- voltage, the eigen-energies for the CQD 
become voltage-dependent, AQ(_a) = Sd ± y/ff^ + rfV'^. 
This dependence remarkably modifies the electronic oc- 
cupation behavior in comparison with no shifting case. 
Electrons have more opportunity to occupy the first QD 
for large bias-voltages, leading to decreasing of the oc- 
cupation number of the second QD, /522, with increasing 
bias- voltage as shown in Fig. [4ja), which is reason of the 
appearance of NDC, because the nonequilibrium current 
is dominatively dependent on the occupation number of 
the second QD in the series configuration. 
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IV. CQD IN A PARALLEL CONFIGURATION 

Focusing attention on the case of resonant tunneling 
through a paraUel CQD, we again choose the optimal res- 
onant condition, 6 = 0, and we also assume the tunneling 
constants as Tlh = Tr22 = T and Tl22 = ^mi = T', 
thus ri^i2 = Tri2 = Vrr', which yield 



^ iix 

^ 22x 



12x 



21x 



r7|(Ax), 



(18a) 

r'/±(AJ+r/±(AJ, (18b) 
rF[e^'^/'/f (Ax) + e-'^/VI (Ax)], (18c) 
rF[e-^'^/Vf (Ax) + e^'^/'/l (Ax)]. (18d) 



At first, we consider the case of ?7 ^ oo. If dot-dot 
hopping vanishes, i.e., fl = 0, the QREs coincide with 
our previous results in Ref. l46i 



A. Simplified QREs in the SR for configurations 
(c) and (f) 

The energetic configurations for coherent tunneling 
through a CQD interferometer are also depicted in Fig. [31 
Analogous to the series CQD, only configurations (c) and 
(f ) bear nonzero tunneling current for the present paral- 
lel arrangement. For reference purposes, we therefore 
provide the simplified QREs written in the SR at zero 
temperature for these two configurations below. 

For an extremely large bias- voltage [configuration (c)], 
these tunneling rates simplify further as ^ii^f^^-^ — 



y- — r r+ 

22q(/3) — ' 22q(/3) 



/pp7gi¥'/2 Y 



12a{f3) 



r+ 

^ 21a(/3) 



^ ' 12q(^)) ~ ^ 21q(/3) ~ 

VTre-"^/'^. The sim- 



lla(/3) 

plified QREs in terms of the SR become finally 
Poo 



Pii = 



P22 



Pl2 



T'pii + Tp22 - (F + F')poo 
+VfF(e^'^/2pi2 + e 

Tpoo - F'pii + ifi(p2i 



1 

rVoo - rp22 - ■i^{p2i 

-iVfF(e*'^/2pi2 



P2lj, 
P12) 



rr'(e^'^/Vi2 + e-^'^/V2i) 



P12) 

e-^'^/V2i) 



(19a) 
(19b) 

(19c) 



i2Spi2 — ifl{pii 

+Vfre'^^^poo - 



P22) - -(r-i-r')pi2 

i\/rFe-*'^/2(pii+p22), 



(19d) 



along with the normalization relation, P00 + P11+P22 = 1, 
and P21 = Pi2- expected, we note that the resulting 
QREs, Eqs. (|19ap - (jl9dp . coincide with those of our pre- 
vious deri;ati^or~same system analyzed in Ref. M 
The current becomes 



II = (r + r')poo, 
Ir = -[r'pii + rp22 



rr'(e^'^/V2 + e- 



(20) 

"'^/'P2l)]. 



Moreover, the system of configuration (f) yields 
Pii 



(Fsg - i\/TT'sin6'cos '^)pqo + i^{p2i - P12) 
-{V + Tel + ^\/rF sine cos |)pn 



-\^/f^{e'^'^p^2 + e-'^'^P2l) 
-iVrF52(e-^^/Vi2 + e'^/V2i) 



Fsin0(pi2 -I- P21) 



1 



(21a) 



P22 



(r Cg - -VFr'sin6lcos —)pm - i^{p2i - P12) 



1 



(F + T's'g + - VFF' sin cos ^)p22 



1 



--VFF'(e''^/V2 + e 
-ix/rFc2(e-*'^/2pi2-He^'^/V2i) 
-iF'sin6'(pi2 + P21), 



2- 

iip 12 \ 
^' P21) 



(21b) 



P12 = 2iSpi2 - ii^ipii - P22) 



sm Wjpoo 



i[x/rFsin^?e''^/2 + T'{sl + 1) + T{cl + l)]pi2 

-F[ivrFe*'^/2- i(F + r')si 

_iVrFe-'^/2(pn+P22) 



'VT'sie 



sm Opii 



P22 ~ -Fsin6lp22- 



(21c) 



The simplified current expressed in the ER is given by 

II = 2[c^^' + s^F-^/^Fsin^?cos((^/2)]poo 

-2[c^r + sIt' + \/lT'sine'cos(^/2)]p„„ 

-[^(r - F') sin 61 + slVrre"^^^ - clVfre-"^^^]po,0 

-[i(F - F') sine - clVfTe^'^/^ + VfPe-^'^/^^a, 

(22) 



which is useful for the analysis of the AB interference 
effect in this configuration. 



B. Calculations and discussion: AB interference 
effect between two pathways 

Below, we present numerical analyses for first-order 
resonant tunneling through a CQD interferometer based 
on the QREs of Eqs. ([7a|) - ([7e)) . and the current formula, 
Eq. p2)) . To start, we emphasize the interference ef- 
fect on tunneling due to the additional pathway. As an 
illustration, we carry out calculations by changing the 
tunnel-coupling strength F' of the additional path from 
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(the series configuration) to F (the completely symmet- 
rical parallel configuration), as shown in Fig. [5] for two 
cases, Ed/F = 1.0 and —1.0, respectively, in the absence 
of magnetic-flux, (p = 0. 



E^/r=-i.o 

£2/r=0.5 




FIG. 5: (a,c) the average current, and (b,d) the differential 
conductance vs. bias- voltage for several increasing tunnel- 
coupling strength values F'/F = 0, 0.1, 0.2, 0.5, and 1.0 of 
the additional path without magnetic-flux, ip = 0, for the 
systems with Ed/F =1.0 (a,b) and Ed/F = —1.0 (c,d). Other 
parameters are fi/F = 0.5 and T/F — 0.05. 



We note that the interference effect significantly 
changes transport properties when the system is in con- 
figuration (f) for Ed/F > 0, while it has little influence 
on tunneling when Ed/F < 0. This asymmetry of the in- 
terference pattern is also due to strong interaction. As 
displayed in Fig. [Sl^c), the eigen-level (3 is always singly- 
occupied in configurations (a) and (e), which blocks entry 
of an additional electron to the CQD, leading to vanish- 
ing current, irrespective of whether the additional tunnel- 
path is available or not. Increasing the tunnel rate F' of 
the additional pathway suppresses the current until it 
entirely vanishes in the case of a completely symmetri- 
cal parallel geometry, F' = F, for the bias- voltage win- 
dow 3 > \V/r\ > 1 of the system with Sd/T = 1.0 and 
fi/F = 1/2 [configuration (f)] as depicted in Fig. [^a). 
Accordingly, the differential conductance finally develops 
a two-peak structure located at the bias-voltage values 
separating configurations (f) and (c) with an enhanced 
peak height [Fig. [5Kb)]. This behavior can be intuitively 
interpreted as a result of perfectly destructive quantum 
interference between the additional pathway and the orig- 
inal one for configuration (f). 



C. Finite-bias-induced AB oscillations and 
magnetic-flux-induced negative differential 
conductance 

In this subsection we examine the magnetic-fiux depen- 
dence of coherent tunneling through the parallel CQD. 
Above, we found that nonzero tunneling current occurs 
only for configurations (f) (finite bias-voltage) and (c) 
(infinite bias- voltage). We consider these two cases. 
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FIG. 6: Calculated magnetic-flux dependence of the dot oc- 
cupations pii and P22 (a) and the tunneling current (b) for 
a CQD with F'/F = 0.5 and Q/F = 0.5 in configurations (c) 
(thin lines) and (f) (thick lines), respectively. 

We determine the nonequilibrium dot occupations 
Pii(22) and tunneling current / based on the simplified 
QREs of Sec. IIV Al The calculated results are illustrated 
in Fig. El for a CQD with F'/F = 0.5 and ^/T = 0.5. 
Clearly, the dot occupations exhibit periodic oscillations 
with period 47r for both cases. Actually, our previous 
studies have already shown that interdot hopping yields 
periodicity of the AB oscillation as AttJ^ We further ob- 
serve that pii ~ sm(p/2) and p22 ^ — sin((^/2) for con- 
figuration (c) , leading to a magnetic-flux independent re- 
sult for Poo = l-2(pii-|-p22)- Since I/T ~ poo [Eq. 
the current of configuration (c) is consequently indepen- 
dent of magnetic-flux, as shown in Fig. EKb). This result 
is quite surprising. In fact, if dissipative mechanisms are 
considered, P22 may not be completely out of phase with 
pii (i.e., partial interference), thus the current of configu- 
ration (c) could become a periodic function of magnetic- 
fiux. 

The current of configuration (f) is complex according 
to Eq. (P^ . The results of the preceding subsection show 
that perfectly destructive quantum interference causes 
the current to vanish for zero magnetic-fiux, </? = 0, in 
this configuration [Fig. [5Ja)]. Application of a magnetic 
field will induce variation of the interference effect from a 
destructive pattern to a constructive one. Therefore, one 
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can expect that the enclosed magnetic-flux, if = ±27r, 
will enhance the current due to perfectly constructive 
interference. A striking result we obtain is that the en- 
hanced current may be much larger than the constant 
current of configuration (c) , indicating the appearance of 
a NDC at the boundary between configurations (f) and 
(c). The variation of the I-V characteristic and differ- 
ential conductance with varying magnetic-flux are shown 
in Fig. m 




FIG. 7; Calculated tunneling current I (a), and differential 
conductance |i (b), for a CQD with F'/F = 0.9, Ed/F = 1.0 
and fi/F = 0.5 for various magnetic-fluxes ((^/ti" = 0, 0.5, 1.0, 
1.5, and 2.0). The temperature is T/F = 0.05. 



D. Finite interdot Coulomb repulsion U 

Here we examine the transport of a CQD in the case 
of weak dot-dot Coulomb interaction, U — 2T. 

Figure Wi^) plots the calculated current for a series- 
CQD with Ed/F = 1.0 and various interdot hoppings, 
n/F = 0.25 ^ 1.0. Different from the case of infinite U , 
a NDC is observed around 2{Xp -I- U) even for a series- 
CQD with weak interdot couphng n/T < 0.5. This NDC 
behavior was first reported in our previous study^^ and 
was ascribed to decoherence effect in coherent tunneling 
due to the coupling to leads. From Fig.[5]^b) and (c), 
we observe that this NDC is robust against the additional 
pathway F' in the absence of magnetic- flux and it is tuned 
to positive differential conductance (PDC) when apply- 
ing the magnetic field. It is worth noting that in the case 
of weak interdot Coulomb repulsion, above mentioned 
magneticc-flux-controUed NDC still remain around the 
boundary between configurations (f) and (c) [here con- 
figuration (c) is at the bias voltage region 3 < V < 5]. 
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FIG. 8: Bias-dependent tunneling current / for a series-CQD 
with various dot-dot hoppings (a), and for a parallel-CQD 
as functions of F' (b) and magnetic-flux (c). The interdot 
Coulomb repulsion is set as U = 2F, and other parameters 
are e^/F = 1.0 and T/F = 0.05. 



V. ZERO-FREQUENCY SHOT NOISE 

A. Two-terminal number-resolved QREs and 
MacDonald's formula 

This section is concerned with the discussion of zero- 
frequency current noise of the CQD interferometer in the 
case of C/ — > oo. For this purpose, we employ MacDon- 
ald's formula for shot noised based on a number-resolved 
version of the QREs describing the number of completed 
tunneling events4^ This can be straightforwardly derived 
from the estabhshed QREs, Eqs. We intro- 

duce the two-terminal number-resolved density matrices 
p^^^\t), representing the probability that the system is 
in the electronic state |x) (for X — x') or in the quantum 
superposition state (for x 7^ x') at time t together with 
n{m) electrons occupying the left (right) lead due to tun- 

nehng events. Obviously, Pxx'(*) = Y.n,mPxx^\^) and 
the resulting two-terminal number-resolved QREs for ar- 
bitrary bias-voltage and interdot hopping are: 

^ ' ■* = ^['^e(-'^iiQ + ^^2/3) + •^e(-'^ii/3 + ^220) 



Poo 



sm0(F+ „ + F+ „ - r+ , - r+ ,)]p(-'") 



1 



(n — l,m) 



'[^^Rlip + '^l^R.2213 ^ 2 '^^^ ^(^fil2/3 + R21p)\P'pl3 ^ 
-- sin 0{V^^-^^ + '^Llip ~ ^L22a ~ ^ L22f3) 
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/ (n-l,m) (n~l,m)^ 
^ \Pa0 + Ppa ) 

sin 6'(r^^^^ + r^ii^ — r^22Q ~ ^^22/3) 



(n — l,m) 



(n.m— 1) 



(23a) 



^ 2 ['^9-'^12/3 



dr. 



21/3 



-sin6'(r^i^ 



22^ 



M I"- 



(n.m) 



The current of lead i] can be evaluated as 

Ctt — t—*oo 

where 



(23d) 



(24) 



(25) 



(n+l.m) 
XPOO ' 

+[c^r+ii„ + s2r+,2„ + 1 sin0(r+i2a + r+2iJ] 

--Sin0(r,,^-r22^)(p^^ +P^a ) 
2^*8^12/3 '-e-L 21/3-'Pa/3 



1 



('„2-p- _ „2-p- X (n,m) 

>>*e^ 21/3 '-e-^ 12p)Ppa ' 



(23b) 



P^/3 "* — [^eTjii^ + (^1^^2213 ~ 9 sin0(rj;^2;3 + ^^21/3)] 



(n+l,m) 
XPOO 



1 



11/3 

(n.rn+l) 
XPoo' 



[s'e^lip + ^6^220 ~ 9 ^™ ^(^12/3 + ^21/3)1^/3/3 



--sin0(rii„ -r22„)(p^/3 ' + p)ja ) 



:('5erj2Q '^9^2x0)^0 



(n,m) 



/3a 



(23c) 



Pa/3 — *^^Pa/3 

+[-sin6'(r^^j^ + rJii^ -rj22Q -rj22/3) 

^2'^ei^L12a + ^^12/3) + 2*e(-'"i21a + ^L21I3)]P0Q^ ' ^ 
+ [- Sm9{r'^j^^ + r^^-^^ — rj^22a ^ ^R220) 
'n^^ei^RUa + ^R123) + n^oi^R21a + ^R21B)]Poa' 



'2^^ fil2a ^ _R12/3-' 2 
" ' ' 22/3 



o['^e(rilQ + r22fj) + '^6(^11,3 + r22a) 



1 



-2 sin0(r+2„ + r+„ - r+2^ - r2+^)]p5") 
+^[c^rr2„ - s^r2-i„ - i sin0(rri„ - r2-2j]pi"r^ 



is the total probability of transferring n(m) electrons 
into the left(right) lead by time t and — n(m) if 
77 = L{R). It is easily verified that the current obtained 
from Eq. ([24]) by means of the number-resolved QREs, 
Eqs. (P5a|) - ((23dl) . is exactly the same as that obtained 
from Eq. (jl2p . The zero- frequency shot noise with re- 
spect to lead rj is similarly defined in terms of p("'™)(i:) 
as well>2^'29did8 



dt 



(26) 



To evaluate Sri(0), we define an auxiliary function 



(27) 



whose equations of motion can be readily deduced em- 
ploying the number-resolved QREs, Eqs (|23ap - (j23dp . 
in matrix form: G^'it) = MrjCit) + Q,^p{t) with 
G^{t) = {Gl^,Gl^,G}p,Gl^,G}J^ and p{t) = 

(poo, Paa,P p£, P a p,Ppa )'^- Mr, and Qr, Can be read easily 
from Eqs. (|23ap ~ (j23d|) . Applying the Laplace transform 
to these equations yields 



G"{s) - {si - M^r'g^pis 



(28) 



where p{s) is readily obtained by applying the Laplace 
transform to its equations of motion with the initial con- 
dition p(0) = Pst [Pst denotes the stationary solution of 
the QREs, Eqs ([7a|) - ([7el )]. Due to the inherent long-time 
stability of the physical system under consideration, all 
real parts of nonzero poles of p{s) and G^{s) are negative 
definite. Consequently, the divergent terms arising in the 
partial fraction expansions of p{s) and G^{s) as s — > 
entirely determine the large-t behavior of the auxiliary 
functions, i.e. the zero-frequency shot noise, Eq. (|26|) . 

It is worth mentioning that (1) our two-terminal 
number-resolved QREs, Eqs. (|23ap - (j23dp . facilitate eval- 
uation of the bias- voltage dependent zero-frequency shot 
noise for arbitrary interdot hopping; (2) our calculations 
yield 5l(0) = ^^(0). 
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B. Results and discussion 

We confined our studies to the case Sd > 0. The Fano 
factor, F = S{0)/2I, is used as the main tool to classify 
current It should be noted that because we use the for- 
mula F = S{0)/2I to calculate it numerically, the Fano 
factor shown in the figures below is physically meaning- 
ful only when current increases above zero to avoid nu- 
merical divergence. Thus, our calculated Fano factor is 
physically meaningful only for configurations (f) and (c). 

To start, we discuss shot noise for a series-connected 
CQD with symmetric geometry. Flu — Tr22 = F. In 
the case of configuration (c), i.e., the extremely large 
bias-voltage limit with zero temperature, we obtain an 
analytical expression for the Fano factor)^ 



SOa;"* -8x^ + 1 
(l-f 12x2)2 ^ 



(29) 



with X = Q/T. For configuration (f), i.e., with appropri- 
ately small bias-voltage, we have 



8192a;'* + 896x2 ^ 

(21 + 128x2)2 • 



(30) 



The results indicate sub-Poissonian shot noise, i.e. F < 
1, in the effective bias- voltage regime for any dot-dot hop- 
ping strength, as shown in Fig. O If the dot-dot hopping 
strength is set moderately weak, il/T < 0.978 (which can 
be achieved by tuning the coupling potential between the 
two dots via applied gate voltage), we have Ff > F^, oth- 




FIG. 9: Calculated Fano factor, F = S{0)/2I, for the series- 
connected CQD with Ed/r = 2.0 and various values of fl/F 
(3/8, 1/2, 1, and 2). The temperature is T/F = 0.05. 

We also examine the effect of the bias- voltage-induced 
shifting of bare levels on the shot noise in Fig. [TOl It is 
shown that this shifting results in (i) a weak enhancement 
of the shot noise, i.e., a super-Poissonian noise to be a 
companion to the NDC when the eigen-level /3 dominates 
in transport, and (ii) a constant Fano factor F = 1/2 at 
the extremely large bias-voltage region. 

Focusing attention on the more general parallel CQD 
geometry, we discuss the quantum interference effect of 
the additional pathway on shot noise. Our calculated 
results in the absence of magnetic-flux are plotted in 
Figs, [m including the shot noise S{0) (normalized to 






v/r 




FIG. 10: (a) Zero-frequency shot noise S{0) (dashed lines) and 
tunneling current I (solid lines), and (b) the Fano factor F, 
for the series-CQD with Sd/T = 1.0, fi/F = 1/2 and shifting 
factors 77 = and 1/4. The temperature is T/F = 0.05. 



^F), the current / (normalized to I^F) for comparison, 
and the Fano factor. We find that the system exhibits 
a huge Fano factor with increasing tunneling rate of the 
additional branch in the case of configuration (c).2^ At 
zero temperature, we arrive at analytical expressions 



Ic 

Fr. 



4x2F(7 + 1) 



(31) 



(7 + 1)2 + 12x2' 

[(8O72 + 3527 + 80)x^ 
+ (-87^ + 1607^ + 3367^ + I6O7 - 8)x2 
+7^ + 107^ + 3I7'' + 447^ + 317^ + IO7 + 1] 
^h-l)-'[h + ir + l2x']-', (32) 



with X = fl/T and 7 ~ F'/F. It is obvious that (1) in the 
case of a series-CQD, Eq. ([32|) reduces exactly to Eq. (|29|) . 
exhibiting sub-Poissonian behavior; and in contrast, (2) 
Fc > 1 with increasing 71^ 





(a) 




f 






E^/r=2.o 










: n/r=i 








s 




















r — y 
/ ^ 


















FIG. 11: (a) Zero-frequency shot noise S{0) and tunneling 
current /, and (b) the Fano factor F, for the case Ed/F = 2.0 
and fi/r = 1 with various values of F'/F in the absence of 
magnetic-flux. The temperature is T/F — 0.1. 

It should be noted that all calculations in the present 
paper are performed under the assumption of full inter- 
ference between the two pathways and infinite inter-dot 
Coulomb repulsion as well. If decoherence is taken into 



13 



account due to some dissipative mechanisms, the degree 
of interference wiU naturally lessen, thus leading to a 
great suppression of the Fano factor. In particular, in 
the case of full noninterference, the Fano factor reduces 
to a constant F — 10/27, exhibiting sub-Poissonian shot 
noise. On the other hand, it has been reported that the 
noise is always sub-Poissonian in the case of no inter-dot 
Coulomb interaction^? 

The phase effect on shot noise in a CQD interferometer 
tuned by magnetic-flux is of special interest for configu- 
rations (c) and (f). We exhibit the magnetic-flux depen- 
dence of current, shot noise, and Fano facor in Fig. [12] 
with r'/r = 0.5 at zero temperature. Due to interdot 
hopping, both the current and shot noise exhibit periodic 
oscillations with period 47r (as well as the Fano factor) 
for configuration (f). The current nearly vanishes around 
zero magnetic-flux, ip = 0, as indicated in Fig. [Hb) in 
Sec. IIV C( while the shot noise is suppressed much more 
due to the effect of perfectly destructive quantum inter- 
ference. In contrast, constructive quantum interference 
enhances the zero-frequency shot noise more than the 
current, giving rise to super-Poissonian noise _F ~ 5/4 
a.t (p = ±277 even for configuration (f). On the other 
hand, the shot noise for configuration (c) displays dif- 
ferent magnetic-flux dependence from that of configura- 
tion (f). Both S{Q) and F are observed to behave as 
~ I sin((^) I . Interestingly, varying magnctic-flux could 
change the shot noise from super-Poissonian, F ~ 10 at 
If = 2nn {n is an integer), to sub-Poissonian, F ~ 0.33 at 
If — 2{n + l)7r. Actually, we have derived an analytical 
expression for noise in configuration (c) and (p — ±7r: 



Fc = [80a;'* -8(72 + 27+1)2:2 +7 
+47^ + 672 + 47 + 1] [(7 + 1) 



^(33) 



indicating sub-Poissonian noise. It is also of interest to 
point out that at these values of magnetic-fiux [tp around 
2(n + l)7r], the quantmn interference effect induces NDC 
even for a very weak tunneling rate of the additional 
pathway, P'/F = 0.1, as indicated in Figs. [6] and [7] in 
Sec. IIV Cl [roughlv ip > 1.27r], but super-Poissonian noise 
is not necessarily an accompaniment of this magnetic- 
flux-tuned NDC. 



VI. CONCLUSIONS 

In summary, we have analyzed first-order transport 
through a CQD AB interferometer with finite dot- 
dot hopping, determining the I-V characteristics, zero- 
frequency current noise, and their magnetic-flux depen- 
dence. To accomplish this, we established generic QREs 
in terms of the ER employing a quantum Langevin equa- 
tion approach in the weak tunneling limit. These QREs 
are valid for arbitrary temperature and bias-voltage, as 
well as arbitrary dot-dot hopping, improving upon our 
previous derivation which was limited by the restriction 
that inter-dot hopping be much weaker than dot-lead 
coupling, i.e., il ^ P. 





FIG. 12; Calculated magnetic-flux dependence of the zero- 
frequency shot noise S(0), tunneling current I, and the Fano 
factor F = 5'(0)/2J for a system with fi/P = 1.0 and P'/F = 
0.5 in configuration (c) [panel (a)], and (f) [panel (b)] at zero 
temperature. 



We have also derived the current and Schottky-type 
shot noise formulae in terms of the RDM elements. Our 
theory proves that the previous scheme for evaluating 
the frequency-independent part of the "classical" intrin- 
sic shot noise (i.e., all tunneling processes providing con- 
tributions to current always yield positive contributions 
to the Schottky-type shot noise) remains valid in the for- 
mulation of quantum-based rate equations in the ER. 

Employing the QREs derived here, we have system- 
atically analyzed coherent resonant tunneling through 
a CQD in a series and a parallel conflguration, respec- 
tively. By discussing variation of the energetic conflgura- 
tion with increasing bias voltage, we have explained the 
asymmetric transport property in the series CQD. We 
have also examined the effect of the bias- voltage-induced 
shifting of bare levels of the CQD and found the ap- 
pearance of a NDC. For the parallel CQD, our numerical 
results have shown that (1) the current of conflguration 
(c) is independent of magnetic-flux due to the combina- 
tion of full interference and the strong Coulomb blockade 
effect; (2) AB oscillations emerge in the current in the 
small bias- voltage regime, i.e., the case of conflguration 
(f); (3) the current nearly vanishes completely around 
ip — due to perfect destructive interference, while it 
is greatly enhanced around ip — 2Tr due to constructive 
interference, and it may even be much larger than the 
current of conflguration (c) , suggesting the possibility of 
magnetic- flux-controllable NDC. 

Finally, we have investigated zero-frequency shot noise 
using MacDonald's formula by rewriting the fully devel- 
oped QREs in terms of a two-terminal number-resolved 
density matrix form. The main result we obtained is 
that the combined effect of interference between two path 
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branches and the infinite inter-dot Coulomb interaction 
may induce a huge Fano factor, which can also be con- 
trolled by manipulating magnetic-flux. 
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